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low-aspect  wing  in  the  limited  flow  of  a nonviscous  fluid 

V.  I.  Kholyavko,  Yu.  F.  Usik 

Aerodynamic  characteristics  are  determined  for  a thin  plane 
wing,  moving  near  a solid  or  a free  surface  (hydrofoil).  In  this 
case  the  theory  of  a thin  body  is  used.  According  to  it  the  spatial 
flow  around  a thin  body  which  is  extended  in  the  direction  of  move- 
ment is  replaced  roughly  by  a two-dimensional  flow  in  transverse 
planes.  In  the  case  under  consideration  the  task  is  reduced  to 
the  study  of  the  movement  of  a plane  plate  near  an  interface 
(.Figure  1) . 

Assume  a plane  low-aspect  wing,  the  maximum  span  of  which 
coincides  with  the  trailing  edge,  is  moving  at  a small  angle  of 
attack  a with  a constant  velocity  Vm.  During  the  time  dt  a sector 
of  the  wing  dx1=V00dt  passes  through  the  fixed  transverse  plane  ZOY. 
In  the  equivalent  two-dimensional  problem  this  movement  of  the 
wing  corresponds  to  the  vertical  shifting  of  a plane  plate,  moving 
with  a constant  velocity  V0=Vooa,  by  a magnitude  VQdt  and  a change 
in  the  width  of  the  plate  from  l(x^)  to 

'<*»)  + dx,  = /(*,)  + ^ Vmdt, 

where  l(x^)  - local  wing  semispan. 

Flow  in  the  plane  ZOY,  caused  by  the  vertical  shifting  of  the 


1. 


plate  and  a change  of  its  width,  leads  to  a change  in  the  associated 
mass  of  the  plate  by  a magnitude,  corresponding  to  the  increase  in 
the  lift  of  the  wing  on  the  length  dx^ , i.e.. 


Integrating  this  expression  for  the  length  of  the  chord 
O^x-^b,  we  obtain  the  lift  of  the  wing 

P 

V==js[dxi  = yi^ni>).  d) 

where  m(b)  - associated  mass  of  the  plate,  determined  in  a section 
of  the  wing  based  on  maximum  span. 

The  magnitude  of  the  transverse  aerodynamic  moment  relative 
to  the  axis  OZ^,  passing  through  the  vertex  of  the  wing,  is  calcu- 
lated by  the  formula 


» » 

= - j *•  s dx‘ = - v-#  a *»dm  * 

p 

x [m(b)b  — j m(xft)d2i]. 

b 

(2) 

The  resisting  force  of  pressure  (force  of  inductive  resistance) 
of  the  wing  taking  into  account  the  drawing  in  effect  on  the  leading 
edges  is  determined  from  the  correlation 

X,  = T Va-  (3) 

We  will  introduce  into  the  consideration,  in  place  of  forces 
and  moments  (1)  - (3),  the  corresponding  coefficients: 


Figure  1 . 

Key:  (1)  Interface. 


where  S - wing  area  in  a plane; 


S = 2j/(x1)dx,i 


q«,  - dynamic  head; 


Then  the  aerodynamic  characteristics  of  a low-aspect  wing  can 
be  presented  in  such  a form: 
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If  the  coefficients  of  lift  and  lateral  moment  are  known,  then 
the  position  of  the  center  of  pressure  of  the  wing  relative  to  its 
vertex  in  shares  of  root  chord  can  be  calculated  using  the  formula 

3. 
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Formulas  (4)  and  (5)  are  obtained  from  the  general  dependences 
of  the  theory  of  a thin  body  and  are  suitable  in  all  cases  when  the 
assumptions  of  this  theory  are  valid.  Let  us  note  two  special 
features  which  follow  from  formulas  (4)  and  (5). 

1.  According  to  (4)  the  coefficients  of  lift  and  moment  are 
a linear  function  of  the  angle  of  attack.  The  results  of  experi- 
mental investigations  confirm  the  linear  nature  of  these  dependences 
for  low  aspect  wings  up  to  a ^ 6°  • 

2.  The  coefficient  of  lift  does  not  depend  on  the  shape  of 

the  wing  in  the  plane  z-^l (x]_)  and  is  determined  only  by  the  cross 
section  of  the  wing  on  a section  of  the  rear  edge.  The  experimental 
data  and  calculations  using  precise  theories  confirm  this  special 
feature  also  for  wings  with  an  aspect  ratio  \ 1.5. 

Thus,  in  spite  of  the  limitedness  of  the  case  X+0,  the  theory 
of  a thin  body  leads  to  qualitatively  correct  results  for  wings  with 
a finite  aspect  ratio.  It  can  be  assumed  that  under  certain  condi- 
tions this  theory  will  also  give  satisfactory  quantitative  results. 

It  ensues  from  (3)  and  (4)  that  the  problem  of  determining  the 
aerodynamic  characteristics  of  a low-aspect  wing  in  a limited  flow 
is  reduced  to  finding  the  associated  mass  of  the  plate  near  the 
interface.  For  calculation  of  the  associated  mass  we  will  use  the 
method  of  singularities  and  we  will  distribute  over  the  length  of 
the  plate  within  the  limits  of  -l^z^l  an  upper  layer  with  a con- 
tinuous intensity  y(z).  It  is  evident  that  in  this  case 
y(z)=-y(-z).  We  will  take  into  account  the  influence  of  the  inter- 
face using  the  method  of  mirror  reflection  (Figure  1). 

The  distribution  of  y(z)  should  satisfy  the  following  condition 
on  the  plate:  the  particle  of  liquid,  adjacent  to  the  plate  in  any 
point  S,  should  possess  a constant  vertical  velocity  Vn  (Figure  1). 
Using  this  boundary  condition,  we  obtain  an  integral  equation  for 
determination  of  the  unknown  function  y(z): 


fund*  f !<*>(*- *.)“* 

),r=T  ± J,  (.-».)•  + 4*1 


2*  Vt. 


(6) 
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Here  the  "plus"  sign  refers  to  the  movement  of  a wing  under  a 
free  surface  (hydrofoil),  and  the  "minus"  sign  - to  the  movement 
of  a wing  near  a solid  surface  (ground) . Values  1 and  h are  determ- 
ined In  the  lateral  plane  of  flow  ZOY  depending  on  the  coordinate 
x^,  which  enters  Into  equation  (6)  as  a parameter.  For  the  calcu- 
lations of  the  coefficients  of  lift  and  inductive  resistance  of  a 
wing  the  values  of  1 and  h should  be  determined  in  the  section  of 
maximum  span. 

We  will  introduce  a new  variable  from  the  correlation 
2=1  cos  (when  -l^z^.1  we  have  ttc-i^cO)  and  make  the  substitution 
y(z)=y(1  cos  ■$‘)=y('$0  , having  designated  h = h/b . Then  equation 
(6)  is  written  as: 


Jcos*-cos»,  J (COS#_co.*,)».ni7  “ " (7) 


Let  us  note  the  partial  solution  of  equations  (6)  and  (7), 
which  corresponds  to  movement  of  a wing  in  a limitless  fluid  (h=<=°). 
From  (6)  and  (7)  with  h=“  we  obtain 


from  which 
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CO*  t — CO* 


db 


2icV„ 


T = 2V.7=L=  = 2V',ctg*.  ^ 

For  the  solution  of  equation  (7)  in  a general  form  with  h/°° 
we  present  the  unknown  function  y(^)  as  a series 

7(#)  = 2V.[^.ct6»+X  ^sin2n*].  (9) 

in  which  the  first  term  when  AQ=1  is  determined  by  the  solution  of 
(8),  and  the  second  term  and  AQ/1  characterize  the  additional  load, 
developing  on  the  plate  from  the  influence  of  the  flow  boundaries. 


It  is  evident  that  when  h-*-°°  the  coefficient  Aq+1,  and  A2n->-0  (n=l, 

2.  . . ) . 

In  the  writing  of  series  (9)  the  condition  y (^)  =-y (it—  v’u)  was 

used . 

If  the  solution  of  (9)  is  known,  then  the  associated  mass  is 
calculated  in  the  following  manner.  With  an  accuracy  to  constant, 
which  subsequently  is  not  essential,  the  values  of  the  function  of 
the  velocity  potentials  are  determined  on  the  surface  of  the  plate 
when  Vq=1: 


<p».  ■ — ±-y  $ 1 (*)  — T -j  J T (®) 


where  the  upper  sign  refers  to  the  upper  surface  ),  and  the 
lower  - to  the  lower  Cf H ) • The  general  formula  for  determination 
of  the  associated  mass  has  the  form  [1] 

e 

Since  in  this  case  on  the  lower  surface  of  the  plate  2-  + >. 
and  on  the  upper  ^ = — l,  and  furthermore  dS=dz  , then 

* 

r 

m=m  —4P j ?■(*>*  = — 4p / J (p,,(*)sln*<**=i 

V 

T 

= 2 p/*  jsindd»fT(*)Sin»dft. 

Substituting  into  this  formula  the  expression  y(£  ) from  (9), 


we  obtain 


m — *p  l*  \A%  + 


Thus  for  determination  of  the  associated  mass  of  the  plate  it 
is  necessary  to  know  the  first  two  coefficients  AQ  and  A 2 of  ex- 
pansion (9).  The  remaining  coefficients  A2n  (n=2,  3,  4...) 
characterize  the  distribution  of  load  on  the  span  of  the  plate  and 
do  not  have  a direct  influence  on  the  total  aerodynamic  character- 
istics of  the  wing. 


6. 


' -1'"  r" 


When  h=“  from  the  solution  of  (8)  it  follows  that  AQ=1  and 
A2n=0,  therefore  according  to  (10) 

m.  = up/*.  ( 11 ) 

The  aerodynamic  characteristics  of  a wing  with  a calculation 
of  (11)  are  determined  using  formulas  (4)  and  (5) 


^ _ w (t>\ 

Here  s - aspect  ratio  of  the  wing 

z^=l(x-^)  - equation  for  the  form  of  the  wing  in  a plane. 
Formulas  (12)  are  the  known  correlations  of  a low-aspect  wing 
in  a limitless  fluid  [2]. 

We  will  switch  to  the  calculation  of  coefficients  of  expan- 
sion A2n  (n=0,  1,  2...).  Equation  (7)  with  a calculation  of  (9) 
is  transformed  as: 

•o 

Ao  — 1^2*  cos  2nd.  ± /(h,  ft,)  =,  1,  (13) 


where  the  following  designations  are  introduced 

i(h,  fts)  = ^(j/,  — cos  ft../,  + ly,)  — 

— y COS  ft,  V]  Ain  (Jiit-i  +-/2i»+l)  + -T  V A 2„  ( J — J tn+i)\ 

n—l  <i—l 

« 

j J_  j" cos  m > d ft 

" * J (cos  d — cos  d“)*  -J-  4fi’  ’ 

After  the  calculation  of  integrals  J and  further  conversions 

m 

of  formula  (13)  we  finally  obtain  the  equations  for  the  coefficients 
of  expansion  A2n  (n=0,  1,  2,...). 

1.  For  movement  of  a wing  near  a solid  surface 

ft,)  + 1 = 2 <4j«  [C2<i(A,  ft.)  — cos  2n  ft,]. 

7. 


(14) 


2.  For  movement  of  a wing  under  a free  surface  (hydrofoil) 


[2  + C,  (*.  ft.)]  - 1 - E A*n  [c„  (A.  ».)  — cos  2/i  ft,]. 


(15) 


The  following  designations  are  introduced  in  (14)  and  (15) 

».)  = — p-|cosd,|  — 2fiy; 

« = — (sin*  0, -f- 4/T*)  ; 

« = y ^ V{+  sin*  ft*  + 4A* ; 

V = V (sin’tt,  + 4 A*)*  + 16  A*  cos*  t>.; 


t = 


77-  “ _ ; T3nlt)  = cos  (2n  arc  cos  /), 

V u*  + 4/1*  • ' 


where  T2n(t)  ~ Chebyshev  polynomial  of  I family. 

If  in  correlations  (14)  or  (15)  a series  of  values  of  is 
fixed,  then  we  will  obtain  a system  of  algebraic  equations  for  the 
determination  of  the  coefficients  of  expansion  A„  (n=0,  1,...). 

In  this  case  the  number  of  fixed  points  . (i=l,  2...)  de- 

si 

termines  the  number  of  equations  in  systems  (14)  and  (15)  and, 
consequently,  the  number  of  unknown  coefficients  A„  (n=0,  1,  2... 
...1-1).  The  remaining  coefficients  (n=i,  i+1,...)  should  be 
accepted  as  equal  to  zero. 

Using  (14)  and  (15)  calculations  were  made  for  the  coefficients 
in  the  case  of  fixing  of  one  point  (I  approximation),  three  (II  ap- 
proximation) and  five  (III  approximation). 

Part  of  the  calculations  for  movement  of  a wing  near  a solid 
surface  are  presented  in  Figure  2.  Based  on  the  known  coefficients 
Aq  and  A^  and  formula  (10)  it  is  possible  to  find  the  associated 
mass  of  a plate,  and  with  formulas  (4)  and  (5)  - the  aerodynamic 
characteristics  of  the  wing. 


Figure  2 . 

Key:  (1)  -approximation;  (2)  Solid  surface. 


Figure  3 gives  the  magnitudes  of  the  derivative  of  the  coef- 
ficient of  lift  based  on  the  angle  of  attack  , related  to  the 
value  of  (?  y a*  in  limitless  flow  (12).  Here  also  the  dots  designate 
the  average  values  of  C £ , obtained  using  the  theory  of  a vortical 

carrying  surface  for  rectangular  wings  with  X=0.4-1.6  [4]. 

Analysis  of  the  calculations  made  shows  that  satisfactory 
results  right  up  to~h^  0.1  can  be  obtained  already  in  the  III  approx- 
imation  for  a wing  close  to  a solid  surface,  and  in  II  - for  a jj 

hydrofoil . 

0 

If  we  are  limited  by  the  values  h } 0 . 4 (wing  close  to  the 

ground)  and  h^,0.25  (hydrofoil),  then  the  aerodynamic  characteristics 
of  a low-aspect  wing  can  be  determined  using  the  I approximation. 

In  this  case  we  have  simple  analytical  dependences,  which  follow 
from  (14)  and  (15)  when  = 7T/2  and  A =0  (n=l,  2,  3...). 
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1.  For  a wing  close  to  a solic  surface 


V 1 + 4i» 


, m = *p  /* 


V l 4-4A* 


Since 


then 


J!L 


Vj+j 4A» 


For  a hydrofoil 


; m = *p/*  ■ 


V\  + 4h' 


Vl  + fr 


V 1 + 4/»* 


Calculations  using  formulas  (16)  and  (17)  are  shown  by  the 
solid  lines. 

V I 4 (2h\% 

Let  us  note  that  if  in  expression  — ; we  replace  2 by  a 

2/i 

larger  number  (3  for  example),  then  formulas  (16)  and  (17)  give 
satisfactory  coincidence  with  calculations  using  the  III  approxi- 
mation up  to  h^O.10. 

Thus  it  is  possible  to  recommend  the  following  approximation 
formulas  for  calculation  of  the  associated  mass  of  a plate  in  a 
limited  flow  when  h^O.10. 

1.  Near  a solid  surface 

V ! + 9^ 

//J=-irp 5 /*. 


2 . Under  a free  surface 


Figure  3- 

Key:  (1)  -approximation;  (2)  Solid  surface;  (3)  Free  surface; 

(4)  according  to  [4]. 


Formulas  (18)  and  (19)  are  suitable  for  the  calculations  of 
moment  characteristics  and  the  position  of  the  center  of  pressure 
of  wings.  For  example,  for  a wing  close  to  a solid  surface  from 
(5)  and  (18)  we  find 


1 — 


bl*(b) 


Before  the  Integral  h=h/l(b),  h - position  of  the  rear  edge  of 
the  wing  relative  to  the  interface,  1(b)  - wing  semispan. 

In  the  integrand  expression  the  magnitudes  of  h and  1 depend 
on  the  coordinate  x^ . At  small  angles  of  attack 

h (*,)  - h + (fc  - x,)  a = / (b)  [S  + ^ (1  - *,)  •] . 

b I 

In  the  assumptions  of  the  theory  of  a thin  body  yjgy  — j-  and 

the  second  term  in  the  right  side  under  certain  conditions  can  be 
of  the  order  of  the  first  term,  therefore  the  position  of  the  center 
of  pressure  on  a plane  wing  close  to  an  interface  depends  on  the 
angle  of  attack. 
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We  will  be  limited  to  the  case  a-*-0,  when  h(x^)=l(b)h=const 
and  the  position  of  the  center  of  pressure  is  determined  by  the 
formula 

Xd^]~ J '* (JC,) V 7* (**>  + 9/? dx'’  (20) 

o 

— I (x  ) 

where  / (x,)  = y~-  - relative  local  wing  semispan;  i,  ^ 

Let  us  consider  a family  of  wings,  in  which  the  shape  in  the 
plane  is  assigned  by  the  equation  /(*,)  — i?- 

The  exponent  m changes  within  the  limits  0^m<i°°.  When  m=l 
we  obtain  a delta  wing,  when  m+O  - rectangular. 

According  to  formula  (20)  we  will  have 

xa  = 1 - 1 gp  [ xinyrxfl  + 9h*dxt.  ( 2 1 } 

V 1 + 91?  ] 


If  similar  transformations  are  fulfilled  for  a hydrofoil, 
then  from  (5)  and  (19)  when  a-»-0  it  follows  that 


xa  = 1 


3 h \ rlfa)  VP(i,)  + 9hl  . 

~y~  FT T*?  J J 2F /*  (*,)  Fs,* — 3*  “ 


(22) 


For  a family  of  wings  /”(*,)  =*xT 


x*  y V »+«*'  J J aVif  + aF-tf  1 


(23) 


When  h=  00  (flow  of  a limitless  fluid)  from  (21)  and  (23) 

».-ETTT-  (2‘0 


This  result  can  also  be  determined  from  formula  (12). 

In  a general  case  the  integrals  in  formulas  (21)  and  (23)  are 
expressed  by  hypergeometrlc  functions  [3].  Figure  ^ gives  the 

calculations  of  *a  =»  ~~  for  wings  with  m=  1/2,  1,  2 (solid  line 
- for  a wing  close  to  a solid  surface,  broken  - for  a hydrofoil). 


Figure  4 . 

An  analysis  of  the  results  obtained  shows  that  with  the  approx- 
imation of  a wing  to  a solid  surface  (ground)  the  center  of  pressure 
is  shifted  to  the  trailing  edge,  and  this  shift  is  greater,  the 
smaller  the  exponent  m.  In  particular,  for  a rectangular  wing 
(m-*0)  the  greatest  shift  should  be  observed. 

The  influence  of  the  boundary  of  a free  surface  for  a hydro- 
foil is  opposite  to  the  effect  of  the  ground. 
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